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ABSTRACT 

Panchromatic observations of the best candidate hyperluminous infrared galaxies from the widest Herschel 
extragalactic imaging survey have led to the discovery of at least four intrinsically luminous z = 2.41 galaxies 
across an ~100 kpc region — a cluster of starbursting protoellipticals. Via subarcsecond interferometric imaging we 
have measured accurate gas and star formation surface densities. The two brightest galaxies span ~3 kpc FWHM 
in submillimeter/radio continuum and CO J = 4-3, and double that in CO J = 1-0. The broad CO line is due 
partly to the multitude of constituent galaxies and partly to large rotational velocities in two counter-rotating gas 
disks — a scenario predicted to lead to the most intense starbursts, which will therefore come in pairs. The disks have 
M dyn of several x 10 11 M©, and gas fractions of ~40%. Velocity dispersions are modest so the disks are unstable, 
potentially on scales commensurate with their radii: these galaxies are undergoing extreme bursts of star formation, 
not confined to their nuclei, at close to the Eddington limit. Their specific star formation rates place them >5 x above 
the main sequence, which supposedly comprises large gas disks like these. Their high star formation efficiencies 
are difficult to reconcile with a simple volumetric star formation law. V-body and dark matter simulations suggest 
that this system is the progenitor of a B(inary)-type «sl0 14 - 6 -M o cluster. 

Key words: galaxies: high-redshift - galaxies: starburst - infrared: galaxies - radio continuum: galaxies - radio 
lines: galaxies - submillimeter: galaxies 
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1. INTRODUCTION 

Of the known denizens of the galaxy zoo, hyperluminous in- 
frared (IR) galaxies (HyLIRGs; L IR ^ 10 13 L G , where the IR 
luminosity is measured across A. rest = 8-1000 /im) are among 
the rarest and most extreme. They provide excellent labora- 
tories with which to confront the most recent hydrodynamic 
simulations of isolated and merging galaxies (e.g., Hayward 
et al. 2011). The Lir of an HyLIRG implies a staggering star 
formation rate, SFR > 1000 M 0 yr -1 (Kennicutt 1998, for a 
Chabrier 2003 initial mass function (IMF)), unless the IMF is 
top-heavy, or there is a substantial contribution to Lir from a 
deeply obscured active galactic nucleus (AGN), in which case 
we are seeing intense star formation accompanied by the rapid 
growth of a massive black hole (Alexander et al. 2008). Either 
way, we are witnessing galaxy formation at its most extreme. 

How best can we identify the most luminous star-forming 
HyLIRGs? A promising method involves searching among the 
brightest sources detected in the widest surveys with Herschel. 
These comprise blazars, low-redshift spirals, and some intrin- 
sically fainter sources that have been strongly lensed (e.g., 
Negrello et al. 2010). However, there is also the intriguing 
possibility that such samples may contain unlensed and thus 
intrinsically luminous galaxies — the rarest HyLIRGs are pre- 
dicted to have a space density, «M0 -7 Mpc -3 and may mark 
the sites of today’s clusters (Negrello et al. 2005; Lapi et al. 
2011 ). 

Of the facilities providing redshifts for bright Herschel 
sources, the 100 m Robert C. Byrd Green Bank Telescope 
(GBT), together with the ultrawide-bandwidth Zpectrometer 
cross-correlation spectrometer and /fa-band receiver, have been 
among the most effective (e.g., Swinbank et al. 2010; Lrayer 
et al. 2011). Having added considerably to the sample of such 
sources with known redshifts, Harris et al. (2012) identified 
a flat trend of L' co with 12 C 16 0 J = 1-0 line 38 width. This 
contrasted with the steep power-law relation between L' co and 
CO line width found by Bothwell et al. (2013) for unlensed 
submillimeter galaxies (SMGs), as expected if the most gas- 
rich galaxies tend to live in the most massive gravitational 
potentials. Harris et al. argued that intrinsically fainter sources 
require a higher lensing magnification, /z, to rise above the 
observational detection threshold for CO J = 1-0, which 
is approximately constant in flux. The flat trend seen in the 
GBT sample is then best interpreted as lensing magnification 
acting on a population with intrinsically steep number counts. 
Those GBT sources closest to the power-law fit seen for 
SMGs — typically those with the widest lines — are likely to be 
suffering the least lensing magnification. This is where we might 
expect to find any HyLIRGs that are lurking among the lensed 
starbursts. 

Our approach here, therefore, is to search for HyLIRGs 
among those bright Herschel lens candidates with the broad- 
est CO lines, starting with the best example in the sample ob- 
served with GBT by Harris et al., HATLAS J084933.4+021443 
(hereafter HATLAS J084933, with 350 /zm flux density, £350 = 
293 mJy), which displays a line consistent with CO J = 1-0 at 
Zlsr = 2.410 ± 0.003, with an FWHM of 1180 ± 320 km s -1 
and Scoi-o = 0-83 ± 0.19 mJy. The amplification predicted 
by Harris et al. for HATLAS J084933 is consistent with unity, 
/z = 2± 1. 


38 Hereafter, CO refers to 12 C 16 0 unless stated otherwise. 


In the next section, we describe an extensive set of obser- 
vations. We present, analyze, interpret, and discuss our re- 
duced images, spectra, and cubes in Section 3, finishing with 
our conclusions in Section 4. We adopt a cosmology with 
H 0 = 71 km s -1 Mpc -1 , = 0.27, and = 0.73, so 
1" equates to 8.25 kpc at z = 2.41. 

2. OBSERVATIONS AND DATA REDUCTION 

Herschel imaging, undertaken as part of the wide-held 
//-ATLAS imaging survey (Eales et al. 2010), led to the se- 
lection of HATLAS J084933 as a potentially distant, lensed 
starburst: distant, because it is an example of a so-called 350 /zm 
peaker, with its thermal dust peak between the 250 and 500 /zm 
bands; lensed, because models of the far-IR/ submillimeter 
(submm) source counts suggest that the majority of objects 
with Ssoo > 100 mJy are expected to be either local (z < 0.1) 
or lensed (e.g., Negrello et al. 2010), along with the occasional 
flat- spectrum radio quasar. 

On the basis that it represented the best chance of finding 
an intrinsically luminous system rather than a lensed galaxy, 
following the arguments laid out in Section 1 , HATLAS J084933 
was selected for observations with the Jansky Very Large Array 
(JVLA). Here, we detail the JVLA observations and those 
obtained with other facilities that followed as a consequence 
of our initial findings, presenting these data sets in order of 
decreasing wavelength. 

2.1. JVLA CO J = 1-0 and 5 GHz Continuum Imaging 

While the National Radio Astronomy Observatory’s 
(NRAO’s) JVLA 39 was in its DnC, C, B, BnA, and A configu- 
rations, between 2012 January and 2013 January, we acquired 
«a30 hr of Ka - band data, scheduled dynamically to ensure ex- 
cellent atmospheric phase and pointing stability. We recorded 
two sets of eight contiguous baseband pairs, 1024 x 2 MHz 
dual-polarization channels in total. We tuned the first set of 
basebands to cover 13 C 16 0 and 12 C ls O J = 1-0. The 12 C 16 0 
J = 1-0 transition (v rest = 115.271203 GHz; Morton &Noreau 
1994) was placed in the second set of baseband pairs, off- 
seting down by 64 MHz to 33.740 GHz to avoid baseband 
edges. 

Around 4 hr of A-configuration C-band data were also 
obtained, recording 1024 x 2 MHz dual-polarization channels 
every 1 s across 4. 2-6. 5 GHz, with a small gap. Typically, 
0.3 GHz was lost to severe radio-frequency interference near 
6.1 GHz, yielding a band center of 5.1 GHz (5.9 cm). 

PKS J0825+0309 and PKS J0839+0319 were observed every 
few minutes to determine accurate complex gain solutions and 
bandpass corrections in the Ka and C bands, respectively. 
3C 286 was observed to set the absolute flux density scale, 
and the pointing accuracy was checked locally every hour. 

The data were reduced using AIPS. The basebands were 
knitted together using the noifs and vbglu tasks, yielding uv 
data sets with two intermediate frequencies, each comprising 
512 x 2 MHz channels. 

Lor the spectral-line data, these were then imaged in groups of 
four channels (71 km s -1 ), with natural weighting. A variety of 
Gaussian tapers were used to weight the data, with distances to 


39 This work is based on observations carried out with the JVLA. The NRAO 
is a facility of the NSF operated under cooperative agreement by Associated 
Universities, Inc. 
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Figure 1. Morphology of HATLAS J084933 at moderate-to-high spatial resolution. Top left: 160, 350, and 870 /xm continuum emission from dust, as detected by 
Herschel/ PACS, CSO/SHARC-2, and SMA, superimposed on a three-color representation of the VISTA z, J, and K s data. We see two bright components: W is 
centered on a red galaxy, near a blue star; T lies at the end of a lenticular galaxy at z = 0.3478, coincident with its red southern tip. Dust emission can also be 
seen coincident with two further galaxies, labeled C and M, to the east of W. Both C and M can also seen in continuum at 5.9 cm (see lower right), at 1.36 mm, 
2.21 mm, and in CO J = 1-0 and J = 4-3. Top right: CO J = 3-2 emission superimposed on a three-color representation of the Spitzer/ IRAC 3.6 and 4.5 /xm 
imaging and a heavily smoothed VISTA J + H + K s image as the blue channel. CO emission is coincident with the two submm-bright clumps, W and T. An apparent 
extension south of T is mi rrored in the IRAC image, suggesting that both are real. Lower left: JVLA CO J = 1-0 imaging superimposed on our HST /WFC3 FI 10W 
snapshot, wherein we see faint diffraction spikes from the stars. Distorted, low surface brightness emission is seen coincident with galaxy C, characteristic of an 
interaction. Lower right: ~ 1"25 x 1"25 stamps centered on M, C, W, and T, showing our JVLA 5.8 cm continuum imaging as blue contours superimposed on the 
HST FI 10W imaging. At the high spatial resolution available here via interferometric imaging, we see a consistent morphological picture for the emission from cool 
gas and dust in the HATLAS J084933 system. The emission is dominated by W and T, which are separated from one another by ~85 kpc in the plane of the sky, 
with significant emission also visible from two fainter galaxies, C and M. Individually, W and T are well resolved spatially. In all panels, contours are plotted at 
—3, 3 x cr, with \/2-spaced increments thereafter, where a is the local noise level. Beam FWHM are shown as solid ellipses. N is up; E is to the left; offsets from 
“2000 = 132?3889 and 8 2000 = 2?2457 are marked in arcseconds. Known stars are labeled. 

(A color version of this figure is available in the online journal.) 


the 30% point of the Gaussian ranging from 100 to 1000 kk, to 
form a number of cubes with FWHM spatial resolutions ranging 
from Of! 53 x O'.' 51 to 2"0 x 1"8. The line-integrated map shown 
in Figure 1 has a \". 1 x 1"0 synthesized beam. 

The naturally weighted 5.1 GHz pseudocontinuum data were 
imaged, averaging all of the channels, loosely following the 
techniques outlined by Owen & Morrison (2008). The resulting 
map has a spatial resolution of 0'. / 53 x (X'42 (P.A., 26°) and an 
rms noise level of 2.9 11 Jy beam -1 (Figure 1). 


2.2. CO J = 3-2 Imaging from CARMA 

We used CARMA 40 to observe the CO J = 3-2 line 
(v rest = 345.795991 GHz, redshifted to v obs = 101.406 GHz). 


40 CARMA construction was derived from the states of Maryland, California 
and Illinois, the James S. McDonnell Foundation, the Gordon and Betty Moore 
Foundation, the Kenneth T. and Eileen L. Norris Foundation, the University of 
Chicago, the Associates of the California Institute of Technology, and the NSF. 
Ongoing CARMA development and operations are supported by the NSF 
under a cooperative agreement, and by the CARMA partner universities. 
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Observations were carried out with the 3 mm receivers and the 
CARMA spectral-line correlator, with an effective bandwidth 
of 3.7 GHz per sideband and a spectral resolution of 15 km s -1 . 
Four tracks were obtained between 2012 February 13 and 17 
in the C configuration (15-352 m baselines), with 14 usable 
antennas and a total on-source observing time of 13.5 hr. The 
nearby quasar, PKS J0757+0956, was used for complex gain 
calibration. The bandpass shape and absolute flux calibration 
(the latter good to ~15%) were derived from observations of 
Mars, PKS J0854+2006 and 3C 84. Using natural weighting, 
we obtained an FWHM beam size of 2"1 x 1"8 (P.A., 126°), 
with a noise level, cr = 0.3 mJy beam -1 (Figure 1). 

2.3. CO J = 4-3 Imaging from IRAM PdB I 

During 2012 February, we obtained data using six 15 m an- 
tennas in the most extended configuration of the Institut de Ra- 
dioastronomie Millimetrique’s Plateau de Bure Interferometer 
(IRAM 41 PdB I), with baselines ranging up to 800 m. 

The observing frequency was set to 135.203 GHz, the 
redshifted frequency of the CO J = 4-3 line (v res t = 
461.04077 GHz). The weather conditions were exceptionally 
good with a precipitable water vapor of around 2 mm and T sys ~ 
100 K on average. PKS J0825+0309 and PKS J0909+0121 were 
used to calibrate the complex gains, while 3C 84 and MWC 349 
were chosen as bandpass and absolute flux calibrators. The rms 
noise level is 0.97 mJy beam -1 in 4 MHz wide spectral chan- 
nels, for a synthesized beam measuring 1"0 x O'/ 5 FWHM (P.A., 
24°). Using only the line-free channels, the rms noise level is 
47 /zJy beam -1 . 

We also report continuum photometry at 1.36 mm from 
observations obtained throughout 2012 as part of a search for 
water that will be reported elsewhere. 

2.4. 870 /zm Continuum Imaging from the SMA 

Approximately 1.4, 4.9, and 1.7 hr of integration were 
obtained using the Submillimeter Array (SMA 42 ) in its compact, 
extended, and very extended array configurations during 2011 
December 8, 2012 January 31, and 2012 March 31, respectively 
(2011B-S044). The receivers were tuned such that the upper 
sideband was centered at 345 GHz (870 /z m), with 8 GHz 
of single-polarization bandwidth in total. 3C 84 was used as 
the bandpass calibrator and Titan was used for absolute flux 
calibration. PKS J0825+0309 and PKS J0909+0121 were again 
used to track and check the phase and gain. 

A robust = 2 weighting scheme resulted in a (X/83 x 
(f/67 beam FWHM (P.A., 86°) and a noise level, o = 
0.77 mJybeam -1 (Figure 1). 

2.5. 350 gim Continuum Imaging with CSO / SHARC-2 

A total of 80 minutes of data were collected at 350 /zm during 
periods of excellent weather between 2012 January 12 and 13 at 
the 10.4 m Caltech Submillimeter Observatory (CSO 43 ), using 
the SHARC-2 camera (Dowell et al. 2003). We took advantage 
of the open-loop, actuated surface of the telescope (Leong 


41 IRAM is supported by INSU/CNRS (France), MPG (Germany), and IGN 
(Spain). 

42 The SMA is a joint project between the Smithsonian Astrophysical 
Observatory and the Academia Sinica Institute of Astronomy and 
Astrophysics, and is funded by the Smithsonian Institution and the Academia 
Sinica. 

43 This material is based upon work at the CSO, which is operated by the 
California Institute of Technology under cooperative agreement with the 
National Science Foundation (AST-0838261). 


et al. 2006) and we checked and updated the focus settings 
several times each night, resulting in near-Gaussian, diffraction- 
limited (8" 5 FWHM) beam profiles. Data were reduced using 
crush (Kovacs 2008) and the resulting maps have a noise level, 
o = 8 mJy beam -1 (Figure 1). 

2.6. Continuum Imaging from Herschel 

The acquisition and reduction of Herschel (Pilbratt et al. 
2010) parallel-mode data from SPIRE (Griffin et al. 2010) and 
PACS (Poglitsch et al. 2010) for the 9 hr Science Demonstration 
Phase field of //-ATLAS (Eales et al. 2010) are described in 
detail by Ibar et al. (2010a), Pascale et al. (2011), and Rigby 
et al. (2011). Here, we use the original SPIRE imaging and 
have obtained much deeper PACS data from the otI program, 
ot1_rivison_1. We recorded data simultaneously at 100 and 
160 /zm while tracking across the target at 20" s -1 in 10 3' 
strips, each offset orthogonally by 4". A total of 180 s on source 
was acquired for each of two near-orthogonal scans. These data 
were tackled with a variant of the pipeline developed by Ibar 
et al. (2010b), reaching cr 4 and 7 mJy at 100 and 160 /zm, 
respectively. The 160 /zm PACS image is shown in Figure 1. 

2.7. Infrared Continuum Imaging from Spitzer, VISTA, 
and the Hubble Space Telescope 

The 3.6 and 4.5 /zm images were acquired using the Infrared 
Array Camera (IRAC; Fazio et al. 2004) on board Spitzer 44 
(Werner et al. 2004) on 2012 June 24 as part of program 80156. 
The imaging involved a 36 position dither pattern, with a total 
exposure time of just over 1 ks, reaching rms depths of 0.3 
and 0.4 /zJy at 3.6 and 4.5 /zm, respectively. Corrected basic 
calibrated data, pre-processed by the Spitzer Science Center, 
were spatially aligned and combined into mosaics (shown in 
Figure 1) with a re-sampled pixel size of 076 and angular 
resolution of 2"-2"5, using version 18.5.0 of mopex (Makovoz 
& Marleau 2005). 

Images were obtained in Z,Y, J,H, and K s as part of 
VIKING, a public survey with the 4 m Visible and Infrared 
Survey Telescope for Astronomy (VISTA) at Paranal, Chile 
(Emerson & Sutherland 2010), with an image quality of (f/9 
FWHM, as shown in Figure 1 . 

A SNAPshot observation was obtained with the Hubble Space 
Telescope 45 (7/ST) on 2012 February 8, as part of Cycle- 19 
proposal 12488, using Wide-Field Camera 3 (WFC3) with its 
wide J filter, F110W. The total exposure time was 252 s. Data 
were reduced using the iraf MultiDrizzle package. Individual 
frames were corrected for distortion, cleaned by cosmic rays and 
other artifacts and median combined. The resulting ~2' x 2' 
image was re-sampled to a finer pixel scale of (X/064 and is 
shown in Figure 1. 

2.8. Optical Spectroscopy from Keck 

We observed several optical/IR sources in the vicinity of 
HATLAS J084933 using the DEIMOS spectrograph on the 10 m 


44 This work is based in part on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology under a contract with NASA. 

45 Based on observations made with the NASA/ESA Hubble Space Telescope, 
obtained at the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS 5-26555. These observations are associated with program 12488. 
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Keck II telescope 46 during 2012 March 2. DEIMOS was used 
with its 600ZD grating (ruled with 600 lines mm -1 ), the GG455 
filter, and a 1" slit at a central wavelength of 720 nm. Two 
positions were targeted, each with a different position angle 
(P.A.), for 40 and 75 minutes. 

3. RESULTS, ANALYSIS, AND DISCUSSION 
3.1. Basic Morphological Description 

At the spatial resolution of SPIRE, HATLAS J08493 is 
essentially unresolved. However, with only modestly improved 
spatial resolution, via deep imaging with PACS at 100-160 /im 
and CSO/SHARC-2 at 350 pun l, the system is resolved into two 
sources, referred to hereafter as W and T (Figure 1). 

Moving to higher spatial resolution, via interferometric imag- 
ing with SMA, IRAM PdBI, and CARMA at rest frame 
255-880 /im, we see a consistent morphological picture for 
the emission from cool dust in the HATLAS J084933 system 
(Figure 1). The dust continuum emission is dominated by W 
and T, which are separated by KX'5, or ~85 kpc in the plane 
of the sky, with W being the brightest of the two from the op- 
tical all the way through to the radio regime. W is among the 
brightest unlensed SMGs ever observed interferometrically at 
submm wavelengths, with T being not far behind (Gear et al. 
2000; Lutz et al. 2001; Younger et al. 2007, 2008, 2009; Wang 
et al. 2007; Aravena et al. 2010; Smolcic et al. 2012a, 2012b). 
Among the unlensed SMG population, its flux density has been 
equaled only by SMM J12371 1+622212 (also known as GN 20, 
Sg -70 = 22.9 ± 2.8 — Pope et al. 2005; Iono et al. 2006). 

Aided by our interferometric astrometry we can see that at 
optical/IR wavelengths the westernmost clump, T, is barely 
visible as a red extension to a lenticular galaxy (Figure 1). W 
coincides with the extremely red component of a red/blue pair, 
separated by ^2". Based on experience with SMGs, such as 
SMM J02399— 0136 (e.g., C. Ferkinhoff et al., in preparation), 
we expected to find evidence of an AGN toward the compact 
blue source, but found instead that it is a star — our optical 
spectroscopy (Section 2.8) reveals absorption due to Mg as well 
as the Na D doublet. 

There is no obvious trend with wavelength for the relative 
IR/submm flux density contributions from W and T (listed in 
Table 1), S W /S T = 1.22±0.33, 2.22±0.42, 1.27±0.19, 1.32± 
0. 17, 1 . 13 ± 0. 10, and 1 .23 ± 0.25 at 100 /mi, 160 /mi, 350 /zm, 
870 /zm, 1.36 mm, and 2 mm, respectively. The error- weighted 
mean for the IR/submm wavelength regime is 1.22 ± 0.07. 
Moving longward to rest frame 1.7 cm, where we expect the 
emission to be dominated by synchrotron radiation related to 
supernova (SN) remnants, HATLAS J084933 is also dominated 
by component W (Figure 1), with 1 S' W /5' X = 2.7±0.7. 

Continuum emission from dust is also seen east of W, from a 
galaxy we label M in Figure 1 . M is detected securely in 870 /zm, 
1.36 mm, 2.21 mm, and 5.9 cm continuum (Table I). 

Still further to the east, emission can be seen from a 
very red, morphologically distorted galaxy — labeled C in 
Figure 1 — which bears a strong resemblance to the Antennae 
(e.g., Whitmore & Schweizer 1995; Klaas et al. 2010). C is 
detected weakly in 870 /zm, 1.36 mm, 2.21 mm, and 5.9 cm 
continuum (Table 1). 


46 Data presented herein were obtained at the W. M. Keck Observatory, which 
is operated as a scientific partnership among the California Institute of 
Technology, the University of California, and NASA. The Observatory was 
made possible by the generous financial support of the W. M. Keck Foundation. 


Table 1 

Continuum Flux Densities 


Wavelength 

5 

W 

V 

T 

Facility 

0.88 /zm 

8.16+1.37 

5.80 ± 0.98 

/zJy; VISTA Z 

1.02 /zm 

10.5 ± 1.8 

8.65 ± 1.45 

/zJy; VISTA Y 

1.1 /zm 

13.6 ± 2.3 

11.2+ 1.9 

/zJy; FI 10W 

1.25 /zm 

16.8 ± 2.8 

13.0 ± 2.2 

/zJy; VISTA J 

1.65 /zm 

23.5 ± 4.0 

19.2 ± 3.2 

/zJy; VISTA H 

2.15 /zm 

38.1 ± 6.4 

26.1 ±4.4 

/zJy; VISTA K s 

3.6 /zm 

67.1 ± 18.6 

47.7 ± 17.4 

/zJy; IRAC 

4.5 /zm 

76.6 + 20.1 

45.8 ± 15.2 

/zJy; IRAC 

100 /zm 

23 ± 7 a 

19 ± 7 a 

mJy; PACS 

160 /zm 

91 ± 12 a 

40 ± 12 a 

mJy; PACS 

250 /zm 

242 ± 18 b 

mJy; SPIRE 

350 /zm 

293 ± 22 b 

mJy; SPIRE 

500 /zm 

231 ± 19 b 

mJy; SPIRE 

870 /zm 

25 + 2 

19 + 2 

mJy; SMA 

1.36 mm 

8.3 ± 0.5 

7.5 ± 0.5 

mJy; PdBI 

2.21 mm 

1.1 ±0.1 

0.8 ±0.1 

mJy; PdBI 

3.0 mm 

0.33 + 0.14 

3er < 0.42 

mJy; CARMA 

8.8 mm 

0.043 ± 0.010 

3a < 0.030 

mJy; JVLA 

5.9 cm 

0.177 + 0.015 

0.065 ± 0.015 

mJy; JVLA 


C 

M 


0.88 /zm 

5.96 ± 1.00 


/zJy; VISTA Z 

1.02 /zm 

5.46 ± 0.92 


/zJy; VISTA Y 

1.1 /zm 

7.66 ± 1.29 


/zJy; FI 10W 

1.25 /zm 

5.72 ± 0.96 


/zJy; VISTA J 

1.65 /zm 

5.09 ± 0.86 


/zJy; VISTA H 

2.15 /zm 

11.1 ± 1.9 


/zJy; VISTA K s 

3.6 /zm 

19.5 ± 5.9 


/zJy; IRAC 

4.5 /zm 

17.7 ± 5.3 


/zJy; IRAC 

870 /zm 

4.6+ 1.9 

6.9 ± 1.6 

mJy; SMA 

1.36 mm 

1.8 ± 0.2 

1.5 ± 0.2 

mJy; PdBI 

2.21 mm 

0.25 ± 0.06 

0.31 ± 0.08 

mJy; PdBI 

5.9 cm 

0.022 ± 0.006 

0.020 ± 0.005 

mJy; JVLA 


Notes. 

a Errors determined by placing many apertures across the map, with 3% and 5% 
calibration uncertainties added in quadrature at 100 and 160 /zm, respectively. 
b Errors include the contribution due to confusion and a 7% calibration 
uncertainty has been added in quadrature (E. Valiante et al., in preparation). 

As we shall see in more detail later (Section 3.5), C and M 
are also detected in our CO J = 1-0 and J = 4-3 spectral-line 
imaging, i.e., they also lie at z = 2.41, alongside W and T. 

Both M and C are coincident with IRAC 3. 6-4. 5 /zm emis- 
sion, suggesting that these are not tidal tails or gaseous 
streams, but rather significant concentrations of stars, gas, and 
dust — luminous SMGs in their own right. 

Our deep rest-frame 133-289 nm Keck spectroscopy 
(Section 2.8), with slits covering W, T, and C, reveals no signif- 
icant line emission from T and C; the only line visible from W, 
albeit very faint, is C ii at 232.6 nm. 

Our CO J = 3-2 image from CARMA is shown in Figure 1. 
It is interesting to note that the low-level southern extension to 
the J = 3-2 emission from T is mirrored at 3. 6-4.5 /zm, which 
suggests that both are real. 

3.2. Leasing Model 

The lenticular galaxy that lies along the line of sight to 
component T lies at z = 0.3478, as revealed by our optical 
spectroscopic observations (Section 2.8), where the slit was 
placed along the major axis of the lenticular. Correcting for 
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Wavelength / m [observed] 

1CT 1 10“ 2 1CT 3 10“ 4 10“ 5 10" 6 



Frequency / Hz [observed] 


Figure 2. Radio-through-optical SEDs of the do mi nant components of 
HATLAS J084933, W (red points, solid line) and T (blue points, dotted line). Our 
measurements cover the SED peak, constraining Lir and 7d us t well. The SEDs 
can be fitted adequately using a combination of thermal dust and synchrotron 
emission components (see Section 3.3) plus the stellar population modeling de- 
scribed in Section 3.4. The key resulting characteristics, Lir, Id us t, Md ust , M stars , 
and <?ir are listed in Table 2. 

(A color version of this figure is available in the online journal.) 

the instrumental resolution, the lines due to Mg and Ca H and 
K suggest the foreground galaxy has a velocity dispersion of 
190 ± 80 km s _1 . Its stellar mass is (1.1 ± 0.4) x 10 11 M Q , 
adopting the approach outlined later in Section 3.4. We have 
modeled this galaxy with an elliptical ( e = 0.3) mass profile, 
truncated smoothly at a radius of 25 kpc (the exact choice has 
little effect since the lensed image is well within this radius). The 
line of sight to the z = 2.41 starburst lies close to the semimajor 
axis of this foreground structure and suffers a magnification 
of /x = 1.5 ± 0.2, where the uncertainty was calculated 
by varying the model parameters — orientation, ellipticity, and 
velocity dispersion — by 10%. 

Galaxy W suffers no significant lensing magnification, as 
far as it is possible to discern from our extensive, wide-held, 
panchromatic imaging data (Section 2.7), which are not easily 
reconciled with a massive foreground galaxy group or cluster. 

3.3. Spectral Energy Distributions of W and T 

Figure 2 shows the spectral energy distribution (SED) of 
HATLAS J084933, concentrating on the brightest components, 
W and T. Our CSO/SHARC-2 350 /rm imaging provides 
evidence that — at the peak of their SEDs — the relative SPIRE 
contributions of components W and T follow the average 
ratio seen in the resolved imaging, 1.22 ± 0.07, as described 
in Section 3.1. Perhaps unsurprisingly, then, we then find 
they contribute in roughly the same ratio, 1.55 dh 0.44, to an 
immense overall bolometric luminosity of HATLAS J084933, 
Lir = (4.8 ± 1.0) x 10 13 L 0 (Table 2), having corrected the 
contribution from T for gravitational amplification (Section 3.2). 

The IR-through-radio photometry can be described ade- 
quately with a model comprising synchrotron and thermal dust 
emission, following Kovacs et al. (2010). Table 2 summa- 
rizes a simultaneous fit to the dominant cold dust tempera- 
ture, L dust , to Lir and q\R (as defined by Helou et al. 1985, 
but where Sir is measured across L rest = 8-1000 /im, as is 
Lir), and to Md us t (for a characteristic photon cross-section to 
mass ratio, = 0.077 m 2 kg -1 — Dunne et al. 2000; Dunne 


& Eales 2001 — where we fixed the frequency dependence of 
the dust emissivity, f, to be 2.0, which is physically plausi- 
ble and consistent with the best-fit value, 2.08 ± 0.15). We 
fixed the synchrotron power-law index, a = —0.75, where 
S v oc v a , and we also fixed y = 7.2 — the power-law index of 
a dust temperature distribution appropriate for local starbursts, 
dMfast/dTfast oc L d ~ s y t , which is designed to offer a physically 
motivated treatment of the Wien side of the thermal emission 
spectrum. Both W and T have 7d ust ~ 36-40 K, commensurate 
with similarly luminous, dusty starbursts at z ~ 2-3, when cal- 
culated using a power-law temperature distribution (Magnelli 
et al. 2012). 

We are sensitive to AGN activity via our optical imaging 
and spectroscopy (Section 2.8) and our radio/mid-IR imaging 
(Sections 2.1 and 2.7). In relation to the far-IR/radio flux ratio, 
W and T have qjR = 2.30 ± 0.08 and 2.53 ± 0.13, respectively, 
spanning (and consistent with) the tight correlation for star- 
forming galaxies (SFGs; q\R ~ 2.4, e.g., Yun et al. 2001). 
Summarizing, we have no evidence 47 that an AGN contributes 
significantly to the Lir of HATLAS J084933, though of course 
we have no evidence of absence since even a powerful AGN can 
be hidden very effectively by the quantities of dust-rich gas in 
this system. 

If contributions to Lir from AGNs are small, then W and 
T are generating stars at a rate of ^3400 and ^1500 M Q yr -1 , 
respectively, for a Chabrier (2003) IMF (Kennicutt 1998, though 
see Dwek et al. 201 1). The estimates of L [R for C and M given in 
Table 2 are scaled from their submm photometry relative to W 
and T, assuming the same SED shape. This suggests that both 
are close to being HyLIRGs in their own right, ^7 x 10 12 L o . 

Before closing our SED discussion, we add a word of caution 
regarding our estimated IR luminosities. If W and T lie in the 
core of a protocluster, as we shall see in what follows, then lower- 
luminosity star-forming cluster galaxies may contribute to the 
flux densities measured by the relatively large SPIRE beam at 
250-500 /xm (see Negrello et al. 2005). Although the calibration 
cannot compare with that of Herschel, SHARC-2 recovers less 
flux than SPIRE at 350 /xm, a possible manifestation of the 
protocluster. 

3.4. Stellar Masses 

The rest-frame ultraviolet-through-near-IR SEDs of W and T 
are illustrated in Figure 2, alongside the far-IR-submm-radio 
data. Flux densities are listed in Table 1 . Our high-quality HST, 
VIKING, and IRAC images at 0.88-4.5 /xm are compromised to 
a certain extent by the unfortunate super-position of a star and a 
lenticular galaxy near W and T, respectively, but the considerable 
uncertainties associated with estimating stellar mass dominate 
over the photometric uncertainties (see, e.g., Hainline et al. 
2011; Michalowski et al. 2012). 

For components W, T, and C, we find best-fitting monochro- 
matic rest-frame H - band luminosities, logL H , of 11.82, 11.46, 
and 11.12 L G , respectively. Stellar masses were determined us- 
ing magphys with the most recent version of the stellar popu- 
lation models of Bruzual & Chariot (2003), a Chabrier (2003) 
IMF, and the SED prior distributions described by da Cunha 
et al. (2008), following Rowlands et al. (2012). The mod- 
els encompass a wide range of exponentially declining star 
formation histories (SFHs), with bursts superimposed (see da 
Cunha et al. 2008 for details). The implied stellar masses are 


47 Observations of higher-7 CO lines would be useful in terms of AGN 
diagnostics. 
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Table 2 

Properties of HATLAS J084933.4+021443 


Property 


W 

T 

M 

c 

R.A. (J2000) 


08:49:33.59 

08:49:32.96 

08:49:33.80 

08:49:33.91 

Decl. (J2000) 


+02:14:44.6 

+02:14:39.7 

+02:14:45.6 

+02:14:45.0 

Mean FWHM at 5.9 cm, 870 /xm, and in CO J 

= 4-3 

2.9 ± 0.4 kpc 

3.8 ± 0.6 kpc a 

b 

b 

FWHM in CO J = 1-0 (kpc) 


7.0 ±2.1 kpc 

6.2 ± 1.4 kpc a 

b 

b 

log L ir (Lq) 


13.52 ± 0.04 

13.16 ±0.05 a 

12.9 ± 0.2 

12.8 ± 0.2 

SFR (M 0 yr- 1 ) (Kennicutt 1998; Chabrier 2003 IMF) 

3400 

1500 a 

800 

640 

Tlust (K) 


39.8 ± 1.0 

36.1 ± 1.1 

b 

b 

<?IR 


2.30 ± 0.08 

2.53 ±0.13 

b 

b 

log M dust (Mq) 


9.32 ± 0.05 

9.10 ± 0.05 a 

b 

b 

CO J = 1-0 7co (Jy km s -1 ) 


0.49 ± 0.06 

0.56 ± 0.07 

0.057 ± 0.013 

0.079 ± 0.014 

CO J = 1-0 L' co (10 9 Kkm s~ J pc 2 ) 


138 ± 17 

157 ± 20 

16.0 ± 3.6 

22.2 ± 3.9 

CO J = 1-0 Leo (10 6 Lq) 


6.72 ± 0.82 

7.65 ± 0.96 

0.779 ±0.178 

0.108 ±0.019 

CO J = 1-0 FWHM (km s -1 ) 


825 ± 115 

610 ±55 

320 ± 70 

250 ± 100 

CO J = 1-0 zlsr 


2.4066 ± 0.0006 

2.4090 ± 0.0003 

2.4176 ± 0.0004 

2.4138 ± 0.0003 

13 C 16 0 and 12 C 18 0 J = 1-0 7 co (Jy km s _1 ) 


3cr < 0.24 c 

3 cr < 0.24 c 

3a < 0.24 c 

3a < 0.24 c 

CO J = 3-2 7 co (Jy km s -1 ) 


4.08 ± 0.92 

5.73 ± 1.13 

0.66 ± 0.27 

1.16 ±0.37 

CO J = 3-2 FWHM (km s- 1 ) 


830 ± 155 

555 ± 75 

260 ± 95 

620 ± 160 

CO J = 3-2 zlsr 


2.4077 ± 0.0009 

2.4096 ± 0.0004 

2.4173 ± 0.0006 

2.414 d 

CO J = 4-3 I co (Jy km s _1 ) 


3.76 ± 0.27 

5.10 ±0.37 

0.66 ±0.10 

0.95 ±0.15 

CO J = 4-3 FWHM (km s -1 ) 


985 ± 115 

545 ± 30 

320 ± 50 

450 ± 90 

CO J = 4-3 zlsr 


2.4068 ± 0.0002 

2.4090 ± 0.0002 

2.4178 ± 0.0003 

2.4149 ± 0.0004 

Lb r 3-2/ i -o 


0.93 ± 0.24 

1.14 ±0.27 

1.29 ± 0.60 

1.63 ± 0.59 

Lb r 4— 3/1— 0 


0.48 ± 0.07 

0.57 ± 0.08 

0.72 ± 0.20 

0.75 ±0.18 

Umax (km s _1 ) 


350 ± 10 

220 ± 15 

20 ± 10 

80 ±20 

R (kpc) 


7.8 ± 1.0 

6.8 ± 0.8 a 

b 

b 

log M dyn (Mq) 


11.51 ±0.05 

11.13 ±0.05 a 

11.11 ±0.15 

11.11 ±0.20 

log A7 H2+ He (Mq), for aco = 0.8 Mq (K km s“ 

-V 2 )- 1 

11.04 ±0.05 

10.92 ± 0.06 a 

10.11 ±0.09 

10.25 ± 0.07 

SFE (Lq Mq 1 ) 


370 

210 a 

760 

440 

log Mstars (Mq) 


11.38 ±0.12 

11.01 ±0.12 a 

~10 

10.36 ±0.18 


Notes. 

a Corrected for /x = 1 .5 ± 0.2. 
b Insufficient S/N to make a useful measurement. 
c Adopting St 0 tai/ Speak as seen f° r 12 C 16 0. 
d Fixed. 


log M stars = 1 1.38+° 0 V 3 , 1 1 -01 + o.n and •0.36+° 0 1 ] 8 7 M© for W, T, 
and C, where the errors are derived from the marginalized prob- 
ability density function, which incorporates the uncertainties 
in the SFH 48 and photometry. Systematic shifts in stellar mass 
of ~0.2 dex result from adopting a Salpeter IMF or ignoring 
thermally pulsating stars on the asymptotic giant branch, and the 
magnitude of photometric contamination by buried AGNs could 
well be similar, so we estimate that our stellar mass estimates 
are accurate to ~0.3 dex. 

Emission from M can be seen at 3. 6-4. 5 [im but the photomet- 
ric uncertainties due to blending with C and W are large so we 
state only that its stellar mass must be considerable, ^10 10 M©. 

There is little sign of the power-law SEDs that can betray the 
presence of AGNs, supporting the lack of AGN spectral features 
reported in Section 3.1, though of course this does not rule out 
the presence of a deeply buried, highly obscured AGN. 

3.5. Spectral-line Characteristics 

Figure 3 shows the JVLA CO J = 1-0 spectra of W and T 
alongside our IRAM CO J = 4-3 spectra. These spectral-line 
data bring us velocity information for the first time, enabling us 
to see that although W and T do lie at approximately the same 


48 Stellar mass estimates are generally robust to changes in SFH (e.g., Ilbert 
et al. 2010, 2013; Pforr et al. 2012). 


redshift and share similar line profiles and 7b ratios, they have 
significantly different line widths (Table 2). Each one is thus a 
distinct, extraordinarily luminous starburst. While the CO line 
profile of T is typical of SMGs (e.g., Bothwell et al. 2013), that 
of W is broader than most, with ~ 1 000 km s -1 FWZI. Both 
profiles are flat-topped or perhaps double-peaked; they cannot 
be described well by a single Gaussian. 

The CO 7 = 1-0 fluxes of W and T measured in our JVLA 
data sum to /coi-o = 1.05 ± 0.09 Jy km s -1 . At first sight 
this is consistent with the Harris et al. GBT measurement, 
1.04 ± 0.37 Jy km s -1 . However, the CO 7 = 1-0 line profiles 
observed by JVLA and GBT appear noticeably different in 
Figure 3: emission in the GBT spectrum continues significantly 
redward (by ~500 km s -1 ) of the emission that can attributed 
to W and T. 

Looking at our IRAM PdBI imaging — see Figure 4 — after 
subtracting the significant continuum emission (several submm 
continuum flux densities are listed in Table 1), we see that 
both C and M are detected securely in CO 7 = 4-3 (at 9.0cr 
and 7.2cr, respectively; we also see 3. 6-4. 5 /rm, submm, and 
radio continuum emission — recall Figure 1). These galaxies are 
therefore unquestionably real and they share approximately the 
same redshift as components W and T. Their positions are listed 
in Table 2, along with those of the two brighter components. 
Their CO emission is centered ^500 km s -1 redward of 
W and T. 
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-2000 -1000 0 1000 2000 

W j=i — o 



-2000 -1000 0 1000 2000 
v /km s -1 (for 2 = 2.4100) 


Figure 3. CO spectra of HATLAS J084933. Upper two panels: CO J = 1-0 
spectra of galaxies W and T superimposed on solid representations of the 
CO J = 4-3 spectra (the former scaled by 16x to put them on the same 
Tj, scale). The two galaxies both have extraordinarily broad CO lines but the 
profiles are noticeably different, so they are not lensed images of the same galaxy. 
The profiles of the CO transitions are indistinguishable for each individual 
component; their Tj, ratios are also similar (Table 2). Middle two panels: 
CO J = 4-3 spectra of galaxies M and C, with line centers significantly redward 
of W and T. Lower panel: JVLA CO J = 1-0 spectra of W and T alongside 
the significantly different profile from GBT/Zpectrometer (Harris et al. 2012). 
Emission at velocities redward of W and T — from M and C, and several other 
faint clumps seen in the JVLA data, spanning a ~100 kpc region — can account 
for the differences between the GBT and JVLA line profiles. 

(A color version of this figure is available in the online journal.) 

As well as C and M, a number of other faint emission features 
can seen in the integrated CO J = 1-0 image (Figure 1). 
Summing C, M, and these faint clumps leads to tolerable 
agreement between the GBT and JVLA line profiles. Although 
we must add 0.87 ± 0.19 Jy km s -1 to 7^, we must also 
increase Iqq T by ~15% to account for the ~8% lower flux 
density assumed for 3C 286 by Harris et al. (2012), and for 
losses due to attenuation of the faint clumps by the GBT primary 
beam. The agreement 49 between the total GBT and JVLA line 


49 Note that the low-level GBT spectrometer baseline (at negative velocities 
relative to HATLAS J084933) may have caused an underestimate of Iqq T . 


intensities remains consistent to within ^1 o and the profile 
measured by the FWHM ^ 22" GBT primary beam is thereby 
reconciled with the JVLA data. 

3.6. Resolved Imaging of the Gas and Dust 

Our interferometric images from JVLA, IRAM PdBI, and 
the SMA have a spatial resolution of ~0'. / 5 FWHM, or ~4 kpc. 
They allow us to compare the sizes of the regions responsible for 
emission from the total molecular gas reservoir (via CO J = 1-0 
emission), from the star-forming gas (via CO J = 4-3), from 
cool dust produced by SNe in regions of recent star formation 
(via submm continuum emission) and from relativistic electrons 
spiraling in or near recent SN remnants (via radio continuum 
emission). 

We find that emission from W and T is resolved in all these 
wavebands, with total flux densities significantly higher than 
the peak flux densities. In Table 2 we quote their deconvolved 
FWHM 50 sizes at 5.9 cm, 870 /xm and in CO J = 4-3 and 
J = 1-0, as determined with two-dimensional (2D) Gaussian 
fits. 

In the absence of significant AGN-powered activity, radio 
continuum observations trace regions where massive stars have 
recently been formed, with no obscuration. Of the tracers 
available to us, we might anticipate the radio emission to be 
the most compact, as indeed is measured. At 870 /xm we see 
no evidence of a more extended emission component, heated by 
older stars — the deconvolved sizes of W and T are consistent, as 
measured in 5.9 cm and 870 /xm continuum, as one might expect 
given the far-IR/radio correlation (Table 2). The CO J = 4-3 
FWHM are slightly larger, as seen for local ultraluminous 
infrared galaxies (ULIRGs; e.g., Wilson et al. 2008), but are 
consistent with the radio and submm values to within the 
uncertainties. 

The error- weighted average of the 5.9 cm and 870 /xm 
continuum and CO J = 4-3 FWHM measurements is (3.4 ± 
0.4) kpc x (2.4 ± 0.4) kpc (P.A., 9 db 15°) for W; for T, before 
correcting for the amplification, we measure (5.6 ± 0.7) kpc x 
(3.5 db 0.7) kpc (P.A., 48 db 16°). In CO J = 1-0, however, W 
and T have considerably larger deconvolved sizes than those 
measured in the submm or radio continuum, as seen for SMGs 
generally (e.g., Ivison et al. 2011; Riechers et al. 2011). 

Summarizing, we have powerful, direct evidence of starbursts 
covering ^3-4 kpc FWHM, with still larger reservoirs of gas 
available to fuel future star formation, on scales of ^6-7 kpc 
FWHM. 

The large sizes measured in CO have been mentioned along- 
side the relatively high L[ Cn ]/L IR ratios measured toward a num- 
ber of SMGs (e.g., Stacey et al. 2010). Both have been taken as 
evidence that SMGs form stars across larger spatial scales than 
the compact, nuclear events seen in local ULIRGs. 

3. 7. Dynamical Characteristics 

The most remarkable characteristics of HATLAS J084933 are 
revealed by the superlative spatial resolution (Section 3.6), and 
velocity resolution and sensitivity of our JVLA CO J = 1-0 
and IRAM PdBI CO J = 4-3 data, as shown in Figure 4. 

To measure the 2D velocity structure, we fit the CO emission 
lines, pixel by pixel, in each of the cubes. Initially, we attempt 


50 We use half-light radii to calculate the area appropriate for gas and star 
formation surface densities, alongside a 0.5 x correction to the SFR or gas 
mass. We use FWHM simply as a convenient way to compare sizes in several 
wavebands. 
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Figure 4. CO J = 1-0 and J = 4-3 imaging of HATLAS J084933. Left: velocity-integrated continuum- subtracted CO emission from our JVLA and PdBI observations. 
In both panels, we mark the positions of the brightest four galaxies, C, M, W, and T. Right: dynamical maps of the galaxies in CO J = 1-0 and J = 4-3. To construct 
these velocity fields, we fit the CO emission in each cube at each pixel, recording the intensity, velocity, and line width where S/N > 5. North is up and east is to the 
left. Inset: one-dimensional rotation curves for W and T, extracted along the major kinematic axis identified in the modeling. Both W and T have rotation curves which 
resemble disks; indeed, W and T are best described as counter-rotating disks, which Mihos & Hemquist (1996) predicted would lead to the most intense starbursts. 
(A color version of this figure is available in the online journal.) 


to identify an emission line in O'/ 2 x 0/2 regions, increasing 
this slowly until we reach S/N > 5, at which point we fit 
the CO profile, allowing the intensity, centroid, and width to 
vary. Uncertainties are calculated by perturbing each parameter 
independently, allowing the remaining parameters to find their 
optimum values, until Ax 2 = 1 is reached. 

We model the 2D velocity fields of W and T by constructing 
and fitting 2D disk models using an arctan function, v(r ) = 
2jt~ 1 u asym arctan(r/r ? ), to describe the kinematics (where 
u a sym is the asymptotic inclination-corrected rotational velocity 
and r t is the effective radius at which the rotation curve turns 
over — see Courteau 1997). We fit six free parameters to the 
disk velocity field: u asym , r t , x, y centroid, P.A., and inclination 
to the sky. The x, y centroid, and r t must be within the field 
of view and the maximum circular velocity must be <2 times 
the maximum velocity seen in the data. Quoted uncertainties 
reflect the range of acceptable models from all of the attempted 
model fits. 

Our velocity cubes reveal that the wide CO line that led to this 
system being selected as a candidate HyLIRG is due partly to 
the multitude of constituent galaxies and partly to large regular, 
systematic rotational velocities in the two brightest constituent 
galaxies. Both W and T are resolved into counter-rotating gas 
disks — the scenario that many studies have predicted should 
lead to the most intense starbursts (e.g., Mihos & Hernquist 
1994, 1996; Taniguchi & Shioya 1998; Borne et al. 2000; Bekki 
2001; Di Matteo et al. 2007; Salome et al. 2012). 


As shown in Figure 5, we find acceptable fits for W and T 
with disk inclinations of 56° ± 10° and 49° ± 10°, at RA.s of 
26° and 145°, with u asym = 415 and 370 km s -1 , respectively. 
Component T is sheared slightly by the foreground lens, as 
evident from the velocity and o maps. Both galaxies show 
small-scale (30-60 km s -1 , rms) deviations from the best- 
fit model. The dynamical masses 51 of W and T are listed in 
Table 2: several x 10 11 M©, where M dyn = R p 2 ax /sin 2 (/). Their 
velocity dispersions, a, corrected for the instrumental velocity 
resolution, are «a65 km s -1 . 

M has no strong velocity gradient. Its u max is 20 km s -1 
(40 km s -1 , peak to peak), o = 140 ± 15 km s -1 for both 
CO J = 1-0 and J = 4-3, and M dyn ^ 1.3 x 10 11 M 0 . For C, 
in CO J = 1-0 we find u max = 80 km s -1 (160 km s -1 , peak 
to peak), cr = 130 ± 20 km s -1 , and M dyn ~ 1.3 x 10 11 M© (in 
CO J = 4-3 we find slightly higher values for a and M dyn ). 

3.8. Star Formation in Disks Off the Main Sequence 

For main-sequence galaxies, star formation is not confined to 
a compact starbursting nucleus (e.g., Elbaz et al. 2011); at z ~ 
2.4, the specific SFR (sSFR) is approximately 2.8 Gyr -1 . W and 
T have sSFRs of 14.2 and 14.7 Gyr -1 , so their starburstiness, 52 
sSFR/sSFR M s > 5, placing them well above the main sequence. 


51 R was taken to be twice the deconvolved half-light radius in CO J = 1-0. 

52 We cannot determine the IRS starburst indicator (Elbaz et al. 2011) for 
HATLAS J084933 since we lack sufficiently deep rest-frame 8 /xm imaging. 
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Figure 5. CO J = 4-3 imaging of galaxies W and T, showing left to right their intensities, resolved velocity profiles, line-of-sight velocity dispersions (er), and our 
best- fit kinematic model (see the text) and the residuals. The color scale for the velocity fields is shown in each panel. We note that there is evidence for lensing shear 
in the velocity field and line-of-sight velocity dispersion for T. 

(A color version of this figure is available in the online journal.) 


As we saw in Sections 3.6 and 3.7 — their star formation is not 
centrally concentrated as would be the case in local starbursts. 
Instead, W and T are undergoing widespread star formation 
in disks that are supported partly by rotation, as envisaged for 
main-sequence galaxies, and share many of the properties of 
massive star-forming disks (e.g., Forster Schreiber et al. 2009). 
Some other factor — presumably their interaction — has taken 
them into the realm of sSFR reserved for starbursts. 

For the Q parameter 53 (Toomre 1964; Goldreich & Lynden- 
Bell 1965), which describes the stability of disks and which has 
been found to tend toward unity in stable situations (in models 
with and without feedback; Hopkins et al. 2012), we find values 
of ^0.35 for both W and T, with a plausible range that falls just 
short of unity. This suggests the disks are unstable and that the 
gas therein is prone to condense and form stars on a timescale 
shorter than the rotational period. 

3.9. Gas Properties and aco 

The brightness temperature (7b) ratios measured for both W 
and T are higher than expected (Table 2), where we usually see 
r 3 _ 2 /i-o ~ 0.52 ± 0.09 and r 4 _ 3 /i_o ~ 0.41 ± 0.07 for SMGs 
(Harris et al. 2010; Ivison et al. 2011; Bothwell et al. 2013), 
which suggests that the thermalized CO J = 3-2 filling factors 
may be close to unity and that we may have found similar sizes 
in CO J = 1-0 and J = 3-2 if our spatial resolution in the 
latter had been two times higher. 

For a CO to H 2 +He conversion factor, aco = 0.8 M© 
(K km s -1 pc 2 ) -1 , as commonly adopted for IR-luminous 
starbursts where the gas is not in virialized individual clouds 
(Bolatto et al. 2013), the gas masses we determine for W, T, C, 
and M are 110, 83.7, 12.8, and 17.8 x 10 9 M© (corrected for 
lensing). 

If we were instead to adopt the all-sample average gas/dust 
mass ratio for nearby galaxies (M gas /Md us t = 9 1 ^ ; Sandstrom 


53 Usually ascribed to a seminal paper on the stability of stellar disks (Toomre 
1964), which contains no mention of Q. 


et al. 2013), we would determine gas masses for W and T of 
19.0 and 1 1.5 x 10 9 M©, yielding gas/dust-motivated aco values 
of 1.4 and 1.1 M© (K km s -1 pc 2 ) -1 , respectively. Due to the 
uncertainty in /cgso, etc., we regard these as accurate to a factor of 
two times at best. For an environment with sub-solar metallicity, 
for example, aco would rise in line with M gas /Md us t- 

A dynamically motivated estimate of aco can be made by 
determining the difference between the dynamical and stellar 
masses in each system, since this represents the total plausible 
mass of gas. 54 Via this method, we find aco = 0-6 and 
0.3 M q (K km s -1 pc 2 ) -1 for W and T, respectively, where the 
uncertainty, ±0.3 dex, is due almost entirely to the notoriously 
difficult estimate of stellar mass (Section 3.4). 

Finally, we can estimate aco using its relationship with gas- 
phase metallicity and CO surface brightness, as developed by 
Narayanan et al. (2012b), exploiting the resolved imaging in 
CO J = 1-0 acquired here. For solar metallicity (Chapman 
et al. 2005) and the sizes that have been measured accurately 
in our JVLA imaging (Section 3.6; Table 2) we find aco = 0.8 
and 0.7 M© (K km s -1 pc 2 ) -1 for W and T, respectively. 

Our mean is aco = 0.8 ± 0.4 M© (K km s -1 pc 2 ) -1 , in line 
with the value usually adopted for starburst galaxies, with a 
range similar to that reported by Downes & Solomon (1998), 
0.3-1. 3 M q (K km s _1 pc 2 ) -1 . To reconcile our dynamical 
measurement of aco with 3.2 M© (K km s _1 pc 2 ) -1 , which 
Genzel et al. (2010) argue is appropriate for all disks, regardless 
of redshift, would require a buried AGN to dominate the near-IR 
light, causing us to overestimate M stars . 

Adopting our best estimate for aco> we find that W and 
T are both approaching the accepted criterion for maximal 
starbursts, with star formation efficiencies, L IR /M gas , of 370 and 
210 L© Mq 1 . C and M are higher still, though with considerable 
uncertainties. The Eddington limit is believed to be around 
500 L 0 Mq 1 ; above this level, the radiation pressure should 


54 Bearing in mind that any AGN contribution to the rest-frame near-IR 
luminosity will led to an overestimate of stellar mass. 
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quickly expel the gas via radiation pressure from young, massive 
stars on dust (Scoville 2004). 

Feedback due to the aforementioned radiation pressure 
will lengthen the total duration of star formation episodes, 
or lead to a series of short bursts. Ignoring this, the gas- 
consumption timescale implied for the galaxies comprising 
HATLAS J084933 is r gas ~M gas /SFR ft* 15-45 Myr, consistent 
with an interaction-driven period of enhanced star formation 
(Scoville et al. 1986), as already hinted by the morphologically 
tortured state of C. 

The stellar masses determined in Section 3.4 for W, T, 
and C suggest gas fractions, / gas = M gas /(M gas + M stars = 
36%-45%. These gas fractions are the top end of the range that 
Narayanan et al. (2012a) showed to be difficult to reconcile with 
cosmological simulations. 

3.10. Relevance for Star-formation Laws 

Data of the quality gathered here have the potential to test 
whether the rules and prescriptions developed to describe star 
formation in the local universe (e.g., Schmidt 1959; Kennicutt 
1989) can be applied to the interstellar medium of distant, gas- 
rich starbursts like HATLAS J084933. Are the data consistent 
with a simple, volumetric star formation law in which the SFR 
is ~1% of the molecular gas mass per local free-fall time, as 
argued by Krumholz et al. (2012)? Are their SFRs more closely 
related to the orbital periods of their entire galactic disks (e.g., 
Genzel et al. 2010; Daddi et al. 2010), or to the quantity of gas 
above some density threshold (e.g., Lada et al. 2010; Heiderman 
et al. 2010)? 

We begin this discussion in the upper panel of Figure 6 : a 
luminosity-luminosity plot of the observables, L IR and L' co , 
for low- and high-redshift SFGs, e.g., z = 0 disks, luminous 
infrared galaxies (LIRGs) and ULIRGs and z = 1-2.5 BzK 
galaxies, taken from Genzel et al. (2010), and for SMGs. For 
this plot, we need make no assumption about the appropriate 
value of oi( 2 Q. The quality and depth of our continuum and 
CO data mean that the errors bars for W and T can barely be 
discerned. Using L' co measurements made only in CO J = 1-0, 
thereby avoiding the need to adopt a T b ratio, alongside self- 
consistent measurements of Lir, we find no compelling evidence 
that different relations apply to SMGs and other SFGs, in line 
with the findings of Ivison et al. (201 1). Here, we use the sample 
of SFGs from Genzel et al. (2010) and find essentially the same 
result. 

It is noteworthy that while component T contributes most to 
Ico or L' co in all three observed CO transitions, component W 
is significantly brighter in the continuum at IR-through-radio 
wavelengths. We find a similar variation between W and T in 
relation to the far-IR/radio flux ratio, with g IR = 2.30 ± 0.08 
and 2.53 ± 0.13, respectively, spanning (and consistent with) 
the tight correlation for SFGs (gi R ~ 2.4, e.g., Yun et al. 2001). 
One might think their different line-to-continuum ratios would 
lead to a relatively large difference in their positions on the plot 
of Lir versus L' co (Figure 6, upper panel). In fact, W and T both 
lie comfortably among the scatter of SMGs, suggesting that this 
scatter may be intrinsic rather than due to measurement error. 

Moving to the standard surface-density Schmidt-Kennicutt 
(S-K) relation, J2sfr oc shown as the middle panel 

of Figure 6, it is now necessary to adopt value(s) for aco 
appropriate for SFGs and SMGs. Genzel et al. (2010) argue that 
3.2 M q (K km s -1 pc 2 ) -1 is appropriate for disks, regardless of 
redshift; our dynamical aco constraints preclude the adoption of 
this value for W and T. We use 3.2 and 0.8 Af© (K km s -1 pc 2 ) -1 




Figure 6. Plots relating the rate or surface density of star formation with the 
available gas. Components W and T are shown as red squares, with error bars, 
in each plot. Circled red symbols denote high-redshift measurements made 
directly in CO J = 1-0. Low- and high-redshift star-forming galaxies (SFGs), 
e.g., z = 0 disks, LIRGs and ULIRGs, and z = 1-2.5 BzK galaxies, are taken 
from Genzel et al. (2010) for all three plots. Top: Lir vs. L' co , where we find no 
compelling evidence that SMGs follow a different relation than the other SFGs. 
Middle: surface-density S-K relation, where W and T lie among hyperluminous 
SMGs, such as SMM J02399— 0136 and GN 20, offset from the relation for 
other SFGs. Lower: Elmegreen-Silk relation between star formation surface 
density and the ratio of gas surface density and dynamical timescale. W and T 
exhibit extreme star formation efficiencies, even relative to SMGs. The fraction 
of their available gas converted into stars per dynamical timescale, e T , is >0.1. 
(A color version of this figure is available in the online journal.) 
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for SFGs and SMGs, respectively, acknowledging that this is 
overly simplistic 55 (e.g., Narayanan et al. 2012b). We find W 
and T nestled among the other SMGs again (Figure 6), but this 
time the SMGs are offset 56 significantly from the relation for 
other SFGs. 

The lower panel of Figure 6 shows the Elmegreen-Silk (E-S) 
relation between star formation surface density and the ratio of 
gas surface density and dynamical timescale. 57 Taken at face 
value, W and T exhibit extreme star formation efficiencies. 
Several other SMGs are equally extreme, converting >10% of 
their available gas into stars per dynamical timescale (e T >0.1). 
Krumholz et al. (2012) have argued in favor of a simple 
volumetric star formation law, linking Galactic clouds, SMGs, 
and all scales between the two, with projection effects leading 
to the observed scatter. Here, we have discovered two disks 
exhibiting star formation of such ferocity that they are capable 
of converting >10% of their molecular gas into stars with each 
rotation. We require a number of parameters to conspire together 
if we are to escape the conclusion that the gas in these disks does 
not obey the simple volumetric star formation law proposed by 
Krumholz et al. (2012). 

Before leaving the S-K and E-S laws, we note that although 
Lir is an effective tracer of XIsfr Eddington-limited star- 
forming disks, Ballantyne et al. (2013) have argued that velocity- 
integrated CO line intensity is a poor proxy for J] gas . Resolved 
observations of the mid-IR rotational lines of H 2 may provide 
the ultimate probe of the star formation relations in high-redshift 
disks, requiring a far-IR space interferometer, such as FIRI 
(Helmich & Ivison 2008). 

3.11. Future Evolution of the System 

What does the HATLAS J084933 system look like in the 
present day? 

The extreme starbursts that brought this galaxy to our atten- 
tion via Herschel are likely to have been triggered by one or 
more interactions (Engel et al. 2010). Certainly, C is distorted, 
morphologically, and the distances between these galaxies are 
consistent with periods of intense star formation during merger 
simulations (e.g., Springel et al. 2005). 

However, without knowledge of the dark matter halo they 
inhabit and the transverse velocities of W, T, C, and M, it is 
impossible to be sure whether these galaxies are gravitationally 
bound, let alone whether they will merge. We find that the 
total energy, 1/2 v 2 + & (km s -1 ) 2 , of each of these systems 
is negative — indicating that the galaxies are bound to one 
another — if their dynamical masses represent as little as 50% of 
the total. This is confirmed by A-body simulations, starting from 
the observed velocities and positions on the sky, and assuming 
that the transverse velocities do not exceed those in the line of 
sight. If these galaxies inhabit an expected dark matter halo of 
mass ^10 13 M q (e.g., Amblard et al. 2011; Hildebrandt et al. 
2013), it is difficult to escape the conclusion that the system is 
bound. 


55 Changing aco by a factor two (in opposite directions) for the SFGs and 
SMGs would bring them in line with one another. Because of the potential for 
a large systematic error in measurements of stellar mass, among other things, 
we view such a change as plausible. 

56 Any change in the area used to calculate the gas or SFR surface density 
results in a point shifting parallel to the power-law relation seen in this plot. 

57 Although we have estimated the orbital timescale, in this calculation, we 
use the dynamical timescale of the disk, calculated using the half-light radius 
and the velocity dispersion, since Krumholz et al. (2012) show that use of the 
orbital time can artificially increase the apparent star formation efficiency. 


Taking another approach, we searched the (500/ h ) 3 Mpc 3 
Millennium Simulation at z = 2 (Springel et al. 2005). In 
each friends-of-friends halo we identified all of the sub-halos, 
as well as the most massive sub-halo, which we take to be 
akin to W. We then computed the velocity and spatial offset 
between the most massive sub-halo and all the other sub-halos 
and tracked the ^200 cases where there are three sub-halos 
within a radius of 300 kpc (comoving) with velocities >90, 
>600, and >800 km s -1 in projection (i.e., taking into account 
likely line-of-sight distances and transverse velocities). Those 
halos, typically «M.6 x 10 13 M 0 at z = 2, end up with a median 
mass of ^3.8+ 5 ’ 5 9 x 10 14 M e in the present day (90% lower limit, 
1.6 x 10 14 M o ) — structures that are thought to have a space 
density of a few xlO -7 Mpc -3 (e.g., Reiprich & Bohringer 
2002; Vikhlinin et al. 2009). 

The scale over which we see the galaxies comprising 
HATLAS J084933, ~ 1 00 kpc, is similar to the core radius of rich 
clusters of galaxies, and the separation of W and T (in terms of 
kpc and km s -1 ) is typical of binary galaxies (e.g., Turner 1976; 
Schweizer 1 987) and within a factor of two times of the proj ected 
separation of the two dominant galaxies in the Coma cluster. 58 
We conclude that such systems can survive for a long time and 
are not necessarily in the process of merging. The four starbursts 
in HATLAS J084933 may instead be signposting a protoclus- 
ter, where the majority of members are likely well below the 
detection thresholds of our panchromatic imaging. It is ironic 
that deliberate attempts to select clusters at submm wavelengths 
(targeting radio galaxies and quasars, for example, Stevens et al. 
2003; De Breuck et al. 2004; Greve et al. 2007; Priddey et al. 
2008) have met with only limited success, while serendipi- 
tous discoveries are becoming relatively commonplace (e.g., 
Riechers et al. 2010). 

3.12. On the Selection of HyLIRGs, and Their Rarity 

Having resolved HATLAS J084933 into its constituent parts, 
if we place W and T individually on the Harris et al. plot of L' co 
versus FWHM line width (Figure 7) we find that W has moved 
further within the scatter of the intrinsic relation, with a predicted 
amplification close to unity. The predicted amplification for T, 
li ~ 2, is consistent with the lens model presented in Section 3.2. 

We can estimate the space density of objects as rare as 
HATLAS J084933 by utilizing the semiempirical model of 
Hopkins et al. (2010) and Hayward et al. (2013). We start with 
halo mass functions and merger rates over the redshift range of 
interest derived from the Millennium Simulation (Springel et al. 
2005), assigning galaxies to each halo following the abundance- 
matching techniques described in Conroy & Wechsler (2009) 
to yield the galaxy-galaxy merger rate as a function of mass, 
merger-mass ratio, and redshift. Utilizing the stellar mass 
functions of Marchesini et al. (2009), with extrapolations to 
higher redshifts following Fontana et al. (2006), we calculate 
the merger rate for all galaxies with a stellar mass ratio equal 
to or greater than that of W and T, with one component at 
least as massive as W. Integrating the merger rates between 
z = 2 and 6 (t = 2.4 Gyr), where we are sensitive to such 
events via A- ATLAS, results in a space density of approximately 
10 7 Mpc -3 , similar to the space density predicted for extreme 
HyLIRGs (Lapi et al. 2011; recall also Section 3.11). For 
the typical gas-consumption timescale found in Section 3.9, 


58 Coma (A1656) is a B(inary)-type cluster, with two bright, supergiant 
galaxies, NGC 4874 and NGC 4889, close to each other in the cluster core 
(Rood & Sastry 1971). 
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Figure 7. L' co vs. CO J = 1-0 FWHM line widths for SMGs with redshifts, 
fluxes, and magnifications from the literature (red squares — for details, see 
Harris et al. 2012, from whence this figure has been adapted). The power-law fit 
to these data represents an estimate of the intrinsic line luminosity vs. line width 
relation for typical SMGs. Blue circles represent the lensed sample from Harris 
et al. (2012). The point for HATLAS J084933 has been replaced with those 
appropriate for components W and T. W has moved well within the scatter of 
the intrinsic relation, with a predicted amplification consistent with unity. The 
predicted amplification for T is consistent with the lens model presented in 
Section 3.2, ji = 1.5 ± 0.2. 

(A color version of this figure is available in the online journal.) 

tsmg ~ 1 5 — 45 Myr, we then expect to find ^3-9 such galaxies 
in the ~100 deg 2 explored so far in //-ATLAS. The vast 
majority of bright Herschel galaxies have been shown to be 
lensed (Negrello et al. 2010), but the low end of this range is 
plausible. This suggests that the most luminous phase of these 
galaxies has a short duration, <15 Myr, and/or that such systems 
require some other rare property — counter-rotating disks being 
one possibility. Stellar mass functions at high redshift have 
significant uncertainties associated with them; nonetheless, this 
calculation emphasizes the rarity of major mergers, such as that 
between W and T, and illustrates the efficiency of the method 
employed here to find extreme HyLIRGs. The observed surface 
density of these systems is such that even with short durations 
for their most luminous phases, they can end up as a subset of 
the massive present-day clusters discussed in Section 3.1 1. 

HATLAS J084933 was selected on the basis of its broad 
CO line profile as a low-magnification, intrinsically luminous 
system, so it is ironic that its wide line profile is in part due 
to the multiple nature of the source and that future searches for 
single maximal-luminosity galaxies may require the selection of 
less extreme profiles. However, fundamentally we conclude that 
the method explored here for the selection of the most extreme 
starburst events in the universe shows promise. 

4. CONCLUSIONS 

1. Exploiting the relationship between CO luminosity and 

line width determined for unlensed starbursts, we have 
identified and removed gravitationally lensed sources from 
the brightest galaxies found in the widest extragalactic 
Herschel survey to yield a sample of intrinsically lumi- 
nous galaxies. Here, we report deep panchromatic follow- 
up observations of the best candidate HyLIRG system, 
HATLAS J084933, with subarcsecond spectral and spatial 
resolution, which have led to the discovery of at least four 


starbursting galaxies across a ~ 1 00 kpc region at z = 2.41, 
each with a significant mass of stars. 

2. The two brightest galaxies, W and T, are separated by 
~85 kpc on the sky. W suffers no gravitational amplifi- 
cation. T is marginally lensed by a foreground lenticular. 

3. The panchromatic SEDs of W and T reveal that both 
are associated with high-luminosity events — HyLIRGs. If 
contributions to Lir from AGNs are small, and to date we 
have found no reason to suspect otherwise, then W and 
T have logL IR ~ 13.52 and 13.16, T dust = 36-40 K, 
and are generating ^3400 and ^1500 M© yr -1 of stars, 
respectively. 

4. Two other gas-rich galaxies, M and C, detected nearby, 
make HATLAS J084933 reminiscent of the first SMG found 
and the best-studied HyLIRG to date, SMM J02399— 0136, 
which comprises a complex galactic nursery with three 
massive stellar components, each seen during a different 
evolutionary stage (C. Lerkinhoff et al., in preparation). 
The four galaxies known to comprise HATLAS J084933 
are most likely bound. A-body simulations and comparison 
with the Millennium Simulation suggest that in the present 
day, the system may resemble a B (inary )-type cluster of 
mass, «M0 14 - 4 * 6 Af 0 . 

5. Subarcsecond interferometric imaging reveals that the two 
brightest galaxies span ~3 kpc LWHM in CO J = 4-3, 
and in submm and radio continuum, and roughly double 
that in CO J = 1-0. Alongside the detection of [C n], this 
supports a scenario involving a widespread burst of intense 
star formation. 

6. Exquisite 3D imaging from JVLA and IRAM PdBI reveals 
counter-rotating gas disks in W and T, a scenario that 
has long been predicted to lead to the most intense bursts 
of star formation. This suggests that similarly luminous 
galaxies will often be found in pairs. Their modest velocity 
dispersions mean that the disks are prone to instabilities 
on many scales, even that of the disks themselves; they 
are undergoing coherent, extreme starbursts, at close to the 
Eddington limit, which are not confined to their nuclei. 
Despite their disk-like morphologies and dynamics, their 
sSLRs place them ~5 times above the main sequence. 

7. W and T have CO line intensities and widths typical 
of the brightest SMGs, with slightly higher 7j, ratios. 
Three independent estimates of the CO-to-H 2 conversion 
factor — exploiting the dynamical mass, the mass of 
dust, and a recent metallicity-dependent relation — suggest 
aco = 0.8±0.4 M q (Kkms -1 pc 2 ) -1 and hence gas masses 
in W and T of ^10 n A/©. Their gas fractions, ~40%, are 
difficult to reconcile with cosmological simulations. 

8. We have determined that the single-dish CO line profile 
of HATLAS J084933 is exceptionally broad because W 
and T have relatively wide lines, broadened further in 
the Harris et al. (2012) integrated GBT spectrum by the 
modest velocity offset between them, and by contributions 
from a number of faint CO clumps slightly redward of 
the two dominant components. We can thus reconcile the 
GBT CO line profile with that of the various galaxies 
found to comprise HATLAS J084933. These situations 
likely happen only rarely, even in the widest extragalactic 
Herschel surveys. We estimate that only a few systems of 
this nature will be found in every 100 deg 2 surveyed. 

9. Placing W and T individually on the Harris et al. plot of 
L' co versus LWHM CO line width, W and T move further 
within the scatter of the intrinsic relation, with a predicted 
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amplification close to unity. We conclude that the method 
outlined here for the selection of extreme starbursts can 
indeed separate HyLIRGs from the more numerous, less 
luminous, strongly lensed population of IR-bright galaxies, 
and may also pinpoint some distant clusters of starbursting 
protoellipticals. 

10. Plotting L ir versus L' co , W and T lie among the SMGs on a 
single relation that is consistent with all SFGs. They also lie 
among the SMGs on the classic S-K and E-S plots, but the 
SMGs now follow a relation that is distinct from those of 
other SFGs. Their star formation efficiency is significantly 
higher, and is difficult to reconcile with a simple volumetric 
star formation law in which ~1% of the available gas is 
converted into stars during each local free-fall time. 

Understanding systems as luminous as HATLAS J084933 
represents a challenge. However, the observations presented 
here have shed considerable light on the characteristics that 
likely lead to the formation — however briefly — of the most 
luminous SFGs in the universe. Observations scheduled for 
Cycle 1 with the Atacama Farge Millimeter/ Submillimeter 
Array promise to reveal more about this dramatic and complex 
protogalactic environment, probing the structure of the disks 
and sensitive to the presence and influence of buried AGNs. 
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